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Abstract 

We study the heat generation in a nano-device with an electric current passing through the device. 
For the first time, a general formula for the heat generation is derived by using the nonequilibrium 
Keldysh Green functions. This formula can be applied in both the linear and nonlinear transport 
regions, for time-dependent systems, and with multi-terminal devices. The formula is also valid 
when the nano-device contains various interactions. As an application of the formula, the heat 
generation of a lead-quantum dot-lead system is investigated. The dc and ac biases are studied in 
detail. We find several interesting behaviors that are unique to nanostructures, revealing significant 
difference from heat generation in macroscopic systems. 

PACS numbers: 65.80.+n, 73.63.-b, 85.65.+h, 71.38.-k 
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In the past two decades, with the development of micro machining technology, the oper- 
ational speed and integrating techniques of semiconducting electronics have made enormous 
progress. At the present time, the size of the electron devices has reached about 90nm in 
mainstream CPU, and its operation frequency has reached several GHz. With these advance- 
ments, one critical issue, the dissipation of the conducting electrons (i.e. the heat generation 
from electric currents), has emerged. This dissipation in nano-devices strongly hinders fur- 
ther development of the semiconductor electronics. Thus, it is important to uncover the laws 
of heat generation induced by electric currents in nano-devices. From a microscopic point 
of view, the main heat generation is due to electron-phonon interactions, through which 
the energy associated with the electric current in an electronic system is transferred to the 
phonon system in the form of heat. The heat generation in macroscopic systems, the Joule 
effect, is well known. However, similar studies in nanostructures are lacking. 

As far as transport is concerned, many studies have been carried out to investigate the 
electronic conduction through nano-devices in the presence of electron-phonon interactions. 
For example, recent experiments have observed the extra conductance peaks in a single-Ceo 
transistoi^ and in suspended carbon nanotubes^ 4 - These extra peaks are due to a phonon 
being either absorbed or emitted while electrons are interacting with the vibrational or 
breathing phonon modes. On the theoretical front, the quantum-dot (QD) systems con- 
taining electron-phonon interactions have been extensively investigated^*^. These studies 
concentrate on physical properties such as the current, the conductance, the noise of the 
current, and the spectral function, etc. Thus, the previous work have mainly focused on the 
electronic transport properties and neglected the heat generation due to electric current. 

In this Letter, we investigate the heat generation while a current flows through a nano- 
device. For the first time, a general exact formula [Eq.(6) or Eq.(8)] of the heat generation is 
derived. It expresses the heat generation in terms of the two-particle Green functions of the 
system, as conventionally done for the conductance formula. This formula is applicable for 
either linear or nonlinear transport, and also for time-dependent systems. The condition for 
zero heat generation is readily achieved with this formula. To illustrate the usefulness of the 
formalism, we investigate the physical properties of heat generation in a single QD system. 
With a dc bias and at zero temperature, a threshold bias V t for non-zero heat generation is 
found. With an ac bias, a finite heat generation is found even at the time when the current 
vanishes. These behaviors are quite different from our knowledge of the heat generation in 
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macroscopic systems. 

Consider the device as a mesoscopic center region coupled to electro-leads that act as 
electron reservoirs (see Fig. la). In the central region, an electron-phonon interaction exists, 
so the heat (i.e. the phonon) is generated when the current flows through the device. This 
system can be described by the following Hamiltonian: 

H = ^u q a\a q + ^\ q {a q A\ + a\A q ) + H e . (1) 
i g 

Here the first term is the Hamiltonian of the free phonon in the central region, with a q (a) q ) 
being the phonon annihilation (creation) operator. The second term expresses the electron- 
phonon interaction and the operator A q = M qii jd\dj, in which d{ is the electron annihi- 
lation operator for the state i in the central region. The last term H e describes the electron 
system and is of the form: 

i a,k 

+ Y^^kiclJi + Hc], (2) 

a,k,i 

where i is the index of the state in the central region, which includes both the energy level and 
spin indices. The center scattering region can be a QD, double QDs, a molecule, a nanotube, 
etc. Moreover, the center region can contain various interactions H int , e.g. the electron- 
electron Coulomb interaction ^ijfery) Uijd\didjdj, the hopping term ^2%Mv^j) Ujd\dj, etc. 
The third and forth terms in Eq.(2) are for the leads and the coupling between the leads 
and the central region, respectively, in which a is the lead index. In Eq.(2), the leads are the 
normal metal or semiconductor. In fact, the derived heat generation formula given below is 
also valid for superconducting and ferromagnetic leads. The device may be under a high- 
frequency ac bias, a time-dependent gate voltage, or a microwave radiation, in those cases, 
the single particle energies e,, e a k, and the hopping element t a ki are all time-dependent.- 

Assume that the isolated phonon system acts as a phonon reservoir and that the phonons 
are in the equilibrium Boson distribution N q = l/[exp(u q /kBT p h) — 1] with the phonon 
temperature T p h& Then the heat generation Q(t) at the time t, i.e. the energy flow from 
the electrons to the phonons, can be calculated from the time-evolution of the operator 
H P h = 52 9 u q a\(t)a q (t): 

Q(t) = (dE ph /dt) = -2Re^w,A,i(aJ(*)^(*)>- (3) 
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We introduce the standard Keldysh nonequilibrium Green functions in the deriva- 
tion: ((X(t)\Y(t>W = -i(T c [X(r)Y(r')]} : ((X(t)\Y(f)})< = -z(Y(t')X(t)), and 
((X(t)\Y(t'))) r = -i0(t - t')([X(t),Y(t')}). Here X and Y represent the Boson opera- 
tors, e.g. a q , A q , djdj, etc., the times r and r' are in the complex-time contour and t, t' are 
the usual time , 7 ' 9 ' 10 and T c is the contour-ordering operator. By using these Green functions, 
the heat generation of Eq.(3) can be expressed as: 

Q(t) = 2ReY,uM(Mt)H(t))) < - (4) 

q 

Since the isolated phonon Hamiltonian has no interaction, one obtains the following equa- 
tions by using the Wick's theorem : 7 ' 9 ' 10 

where g q (r, r') = ((^(t^oJJt')))^ is the free Green function of the decoupled phonon system 
(i.e. when X q = 0). After the Green function on the complex contour is solved, the lesser 
component is simply of the form: 

«Ar(*)|fiI(0» < = \j dh {((AMAlihWg^t) + <<i,(t)|^0>> < S?(*i,O} . (5) 

Here the free phonon Green functions g<(t,t') = —iN q e~ lu)q ^ t ~ t " > and g q (t,t') = i0(—t + 
t')e~ iu3, '( t ~ t '\ Substituting ((A q (t)\a q (t))) < into Eq.(4), the heat generation Q(t) is obtained 
as: 

Q(t) = -2ImJ2^K f cft 1 e^-*){-^<(^(t)|it(t 1 )>>'- + <<^ I (t)|it( tl ))><}. (6) 

q J ~°° 

If the phonon spectrum is continuous and multi-modes, the summation Y2 q i n Eq.(6) changes 
into an integral ^Tg J °° Pp{oJ q ) dui q , with the phonon mode index (3 and the phonon density 
of states pp. In this case, the heat generation Q(t) is: 

Q(t) = -2ImJ2 du qP ^ q \% f ^ 1 e-^^-*){((i^(t)|4 9 (t0)) < -^((i /3g (t)l4 g (ti))) r }(7) 

a JO J-OO k J 

Eq.(6) or Eq.(7) is one of the main results of this paper. Without any approximation, 
it expresses the heat generation in terms of the electronic two-particle Green function for 
the central region, as commonly done for transport coefficients. This formulation can be 
applied to the nonlinear case with a finite bias, to an ac bias case. It is also valid for various 
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interactions possibly existing in the central region, or with superconducting or ferromagnetic 
leads, etc. 

Next, we consider the steady state case when the Hamiltonian (1) is time- independent. 
In this case, the two-particle Green functions in Eq.(6) depend only on the time difference 
t — t\. By introducing the Fourier transform of the Green functions: {{A q (t)\A q (ti))) r,< = 
j ^ e -^(t-ti)^^|^t))r,< ; t hen Eq.(6) reduces to: 

Q = Ys"*K + • ( 8 ) 



If a Hamiltonian H has the property that [A q , H] = 0, then A q (t) = A q (0) and {{A q \A q )) r j< 
for any non-zero u q . Thus, the heat generation is identically zero from Eq.(8). This means 
that the value [A q , H] represents the strength of the dissipation in the central region. Notice 
that although the heat generation is zero for [A q , H] = 0, the electron spectral function as 
well as the electron transport properties are still affected, sometimes even strongly affected, 
by the presence of electron-phonon interactions.^ 

If the interactions in the central region are weak, then the two-particle Green func- 
tions in Eq.(8) can be analyzed by making the approximation (T c [d\(T)dj(r)d k (0)di(0)]) = 
-(T c [c/|(r)^(0)])(r c [4(0)4(r)]).^ Under this approximation, ((i g (r)|it(o))) c = 
- J E ij ,fc,i M 9,#,*iJcG J C f ; (T)G( i (-r) 1 where G?-(r) = ((di{r)\d]{0))) c being the single-electron 
Green function in the central region. Then the lesser and retarded components can 
be easily obtained by using: [G c (t)G c (-t)}< = G<(t)G>(-t) and [G c (t)G c (-t)Y = 
—i[G > (t)G a (—t) + G r (t)G > (—t)} Taking the Fourier transformations and substituting them 
into Eq.(8), we have: 

Q = Re ^K M ^Kikf ^{GfMGfM 

q,i,j,k,l 

- 2N q [G> k (u)Gl(u) + G r jk (u)G>(u)] } , (9) 



where u = uj — uj q . In this expression, the heat generation is described by the single-electron 
Green functions in the central scattering region. In fact, since these single-electron Green 
functions have been obtained in previous transport studies,—^ so the heat generation can 
be calculated straightforwardly by using this expression. 

In order to illustrate power of our formalism, we now study two specific examples. Con- 
sider a lead-QD-lead system where the QD has an electronic level that is coupled to a local 
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single-phonon mode (e.g. the optical phonon 3 , or the vibrational phonon mode^). The 
Hamiltonian of this device is: 



H = uo q a q d q + \ q {d q + a) q )d)d + e d Sd 

+ ^^(t^ikaak + ^UciJ + H.c). (10) 

a,k ce,k 

Here a = L,R represent the left and the right leads. For this Hamiltonian, the operator 
A q = d)d. We consider two specific cases: (i) A dc bias is applied to the left and right leads, 
with the time-independent e a k- (ii) An ac bias (without a dc component) is applied to the 
two leads, with e a k(t) = e° ak ± Acosc<vi. In order to solve for the Green functions, we first 
make a canonical transformation with the unitary operator U = exp{(\ q /u q )(d q — d q )d)d]. 
After this transformation, the Hamiltonian becomes: 

H = u q a q d q + e d d)d 

+ ^e ak (t)al k d ak + ^t a (c ] ak dX + H.c), (11) 

a,k a,k 

where e d = e d — X q /uj q and X = exp{ — (X q /uj q )(d q —d q )}. Under the canonical transformation, 
the operator A q is unchange, i.e. UA q (t)U^ = Ud) (t)d(t)W = d){t)d{t) = A q (t). Next, we 
make an approximation, replacing the operator X in the Hamiltonian (11) by its expectation 
values < X >= exp{ — (\ q /uj q ) 2 (N q + l/2)}. This approximation is valid when t a <C A g . 6ll ° l12 
Under this approximation, the electron-phonon interaction is decoupled, allowing the two- 
particle Green functions in the heat generation formula (6) or (8) to be expressed by the 
single-particle Green functions, thus, the problem is solved in a straightforward way. 

(i) The dc bias case. In this case, the Hamiltonian is time- independent. After the above 
approximation, the heat generation formulation (8) is reduced to: 

Q = ^ J ^{g<mg>( w ) 

- 2iV,Re [G>(u)G a M + GT d (u)G>(u))] } . (12) 

Here Gd are the QD's single-electron Green functions of the Hamiltonian (11), and they 
are easily obtained as: G^ = iG r d [f L f L + f R f R ]G% G^ = -iG r d Y, a ^cx{l - fa)G d , and 
G r d {u) = Gf(u) = l/{u-e d +iT}, where f a (u) = l/{exp[(u - n a )/k B T e ] + 1} is the Fermi 
distribution, T = (T L + T R )/2, and T a = Y a exp[—{X q /uj q ) 2 {2N q + 1)] with the line-width 
functions T a = 2np a \t a \ 2 . Then after substituting these Green functions into Eq.(12), and 



assuming that T a is independent of the energy uj (i.e. the wide linewidth approximation), 7 
the heat generation is obtained as: 



where f a = f a (u> — uj q ) and N e = l/[exp(ou q /kBT e ) — 1]. The heat generation of Eq.(13) is 
separated into two terms. The first term is proportional to N e — N q and it represents the 
heat conduction between the electron system and the phonon system due to the temperature 
difference (T e ^ T p h). If T e > T p h, this term is positive, meaning that the heat flows from the 
electron system to the phonon system. On the other hand, if T e < T p h, this term is negative. 
The second term describes the heat generation due to the electronic current flowing through 
the device. In the case of wide linewidth, this term is always positive regardless of the values 
of fjLji, erf, and the temperatures T e and T ph . In other words, the heat is emitted when 
the current flows through the QD. However, if T a depends on the energy u, this term can be 
negative in some cases. A negative Q value means that the heat is absorbed when a current 
is flowing, making this system work as a refrigerator. 

We now numerically investigate the heat generation Q for symmetric barriers (Tl = 
Fr = T) and with the same electron and phonon temperature, T e = T p h = T. Fig. lb 
shows Q versus the bias V (V = [il — //#) at different temperatures. At zero temperature 
(T = 0), there exists a threshold bias V t = Wq, and the heat generation Q vanishes when 
V < Vt, independent of other parameters. The reason is that the process of emitting a 
phonon requires an occupied state at uj and an empty state at uj — uj q , and this process is 
impossible for V < uj q at T = 0. For non-zero temperatures, the heat generation is still 
quite small so long as V < V t . However, with the bias increasing and passing uj q + 2T, 
the heat generation Q increases rapidly. At a large bias limit, the heat generation has a 
saturated value uj q \ 2 q Y / {uj^ + 4T 2 ). The curve of dQ/dV versus V does not exhibit the 
phonon side-peaks. For a comparison, the current I versus the bias V is also shown in 
Fig.lb.^2 The current has no threshold bias and it is non-zero even for very small bias V at 
T = 0. Moreover, the curve of dl/dV-V exhibits the phonon side-peaks 

Figures lc and Id show the heat generation Q and the current I versus the QD's level 
6jj at T — and for fixed bias V. For V < V t , Q is zero everywhere, but the curve I -id has 




(13) 
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a peak at e d = 0. When V is slightly larger than V t = oo q , the curve Q-e d has two peaks at 
±uj q /2, while the current I is flat for e d within ±.oo q /2. For a large bias V, Q is large when 
e d is around 0. But the shape of the curve Q-e d is different from that of I-e d . This means 
that the heat generation Q is not proportional to I. However, it is well known that in a 
macroscopic system, the Joule heat is proportional to I for a fixed bias V. 

(ii) The ac bias case. When the operator X in the Hamiltonian (11) is replaced by 
< X >, the electron-phonon interaction is decoupled, and the Green functions, then the 
time-dependent heat generation Q(t) and the current I(t), can be easily obtained. Fig. 2 
shows the heat generation Q(t) and the current I(t). The current I(t) is a periodic function 
with the period 2ir/u ac , and the time-averaged current < I(t) > is zero. The time-dependent 
heat generation Q(t) is also periodic, but with a different period 7i/uj ac , and the time-averaged 
heat generation < Q(t) > is always positive. In particular, Q(t) is non-zero even at the time 
t when I{t) = 0. 

In summary, a general formula for the heat generation due to a current flowing in nano- 
devices has been derived. This formula expresses the heat transfer to the phonon system 
in terms of two-particle non-equilibrium Green functions. Employing this formula, we have 
investigated the heat generation in a lead-QD-lead system. With a dc bias, there exists 
a threshold bias V t at zero temperature. For a fixed bias V, the heat generation is not 
proportional to the current /. Under an ac bias, the time-dependent heat generation Q(t) 
is studied and Q(t) is not zero even when the current I(t) is zero. These results are absent 
in macroscopic systems, unique to nano-structures. 
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FIG. 1: (Color online) (a) Schematic diagram for the system of a central scattering region coupled 
to the leads and the phonon reservoir, (b) The heat generation Q (thick curves) and the current 
I (thin solid curve) vs. the bias V at different temperature T with = 0. (c) and (d) are Q and 
I vs Id at T = 0. The other parameters are u> q = X q = 3 and T = 1. 

FIG. 2: The time-dependent heat generation Q(t) (solid curve) and the current / (dotted thick 
curve) vs. the time t with the parameters uo ac = 5, A = 10, = 0.2, uj q = X q = 3, F = 1, and 
T = 0. The thin dotted curve schematically shows the ac-biased voltage 2A cos u; ac i 
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